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ABSTRACT 

Aims. We present a model for dm-fiber bursts that is based on assuming fast sausage magnetoacoustic wave trains that propagate 
along a dense vertical filament or current sheet. 

Methods. Eight groups of dm-fiber bursts that were observed during solar flares were selected and analyzed by the wavelet analysis 
method. To model these fiber bursts we built a semi-empirical model. We also did magnetohydrodynamic simulations of a propagation 
of the magnetoacoustic wave train in a vertical and gravitationally stratified current sheet. 

Results. In the wavelet spectra of the fiber bursts computed at different radio frequencies we found the wavelet tadpoles, whose head 
maxima have the same frequency drift as the drift of fiber bursts. It indicates that the drift of these fiber bursts can be explained by the 
propagating fast sausage magnetoacoustic wave train. Using new semi-empirical and magnetohydrodynamic models with a simple 
radio emission model we generated the artificial radio spectra of the fiber bursts, which are similar to the observed ones. 
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1. Introduction 

Fiber bursts are considered to be a fine structure of the broadband 
type IV radio bursts. In the dynamic spectrum, they occur in 
groups of many single fibers, which are narrowband and have 
almost the same frequency drift. Because this frequency drift is 
between those of types II and III bursts, the fiber bursts are also 
called the intermediate-drift bursts (Bernold & Treumann 1983, 
Aurass et al. 1987, Benz & Mann 1998, Jificka et al. 2001). 

There are two types of models of these bursts: a) one based 
on the whistler waves (Kuijpers 1975, Chernov 1976, Mann et 
al. 1987) and b) one based on the Alfven waves or solitons 
(Treumann et al. 1990). Both these types of models can be used 
for estimating the coronal magnetic fields (Benz & Mann 1998, 
Rausche et al. 2007), which is otherwise very difficult to deter- 
mine. For this reason the theory and interpretation of fiber bursts 
are very important. However, these two types of models give 
different values for the magnetic field strength, and there is still 
a debate about which models are correct. 

Besides these models, Kuznetsov (2006) proposes a model 
in which the fiber bursts are generated by a modulation of the ra- 
dio emission by magnetohydrodynamic waves. He also proposes 
that these waves could be magnetoacoustic waves of a sausage 
mode type that propagates along the dense coronal loop. 

It is known that the structures with higher plasma density 
(like dense loops or current sheets) act as waveguides, where the 
fast magnetoacoustic waves are trapped (Roberts et al. 1984). If 
these waves are impulsively triggered at some location in such 
a dense structure, then at some distance from the initiation site 
along this structure, these propagating waves form the wave 
trains owing to their dispersion properties. The trains exhibit 
three phases: (1) periodic phase (long-period spectral compo- 
nents arrive as the first ones at the observation point), (2) quasi- 
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periodic phase (both long- and short-period spectral components 
arrive and interact), and finally (3) decay phase (Roberts et al. 
1983 and 1984). As shown by Nakariakov et al. (2004), wavelet 
spectra of these trains correspond to the wavelet tadpoles (tail 
plus head). Detailed numerical simulations of these magnetoa- 
coustic (fast sausage) wave trains and corresponding tadpoles in 
the density slab, as well as in the current sheet, were performed 
by Jelinek & Kai'licky (2012). 

Meszarosova et al. (2009a, b) analyzed several radio events. 
In the first paper we found the wavelet tadpoles in the gyro- 
synchrotron radio burst, where all the tadpoles were detected at 
the same time in the whole frequency range. In second paper, 
analyzing the 2005 July 1 1 radio burst (generated by the plasma 
emission processes), we found the wavelet tadpoles that slowly 
drifted with the frequency drift corresponding to the drift of the 
whole group of the fiber bursts. Furthermore, Meszarosova et 
al. (2011) repeated the wavelet analysis of the 2005 July 11 ra- 
dio event, considering a much broader range of time periods, as 
well as longer time intervals than in the previous paper Thus, 
we found the tadpoles with the period of P w 1.9 s and with 
the fast frequency drift corresponding to the drift of individual 
fiber bursts. Based on this, we proposed that the fiber bursts are 
generated by a modulation of the type IV radio emission by the 
magnetoacoustic wave trains. Karlicky et al. (2011) also found 
the wavelet tadpoles in sources of the narrowband dm-spikes. 

In the present paper, we firstly verify the above-mentioned 
statement about the similarity of the frequency drifts of individ- 
ual fiber bursts to head maxima of the corresponding tadpoles 
by the wavelet analysis of more observed groups of fibers. Then 
we simulate fibers in a newly developed semi-empirical model. 
Due to limitations of this model we also made a first attempt to 
simulate the fiber bursts in a full magnetohydrodynamic (MHD) 
model with the fast sausage magnetoacoustic wave train prop- 
agating in a vertical and gravitationally stratified current sheet. 
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Fig.l. Upper panels: Examples of the fiber bursts at 12:01:05-12:01:16UT (November 23, 1998, left panel) and 8:35:05- 
8:35: 17 UT (November 18, 2003, right panel). Bottom part: Corresponding wavelet power spectra showing the tadpoles with the 
period P = 1.4 s (Table 2). In both spectra, at selected frequencies the times of the tadpoles' head maxima were determined (fl-f6) 
and also shown in the upper dynamic radio spectra (upper panels). The tadpole head maxima drift as the fiber bursts. 



Using a simple radio emission model we produce artificial radio 
spectra of the fiber bursts, which are similar to those observed. 

The paper is organized as follows: Section 2 describes the set 
of observed fiber bursts and their wavelet analysis. In Section 3, 
we present our semi-empirical model of fiber bursts and their 
simulations. Section 4 shows a simulation of the fiber bursts us- 
ing the MHD model. Finally, in Section 5 the results are dis- 
cussed. 



2. Fiber bursts observation and its waletet analysis 

We selected eight decimetric radio events with the fiber bursts 
recorded during the years 1998-2005 by the Ondfejov ra- 
diospectrograph (Jificka et al. 1993). The times of their obser- 
vations, the GOES X-ray and Ha characteristics of associated 
flares are presented in Table 1 . Examples of the characteristic ra- 
dio dynamic spectra that show two typical groups of fiber bursts 
are presented in upper panels of Fig. [T](1998 November 23, 
12:01:05-12:01:16UT - left panel, and 2003 November 18, 
8:35:05-8:35: 17 UT - right panel). The parameters of all fiber 
events are summarized in Table 2. The fiber bursts were ob- 
served in the frequency range 1000-2000 MHz, and their fre- 
quency drift FDf ranges from -73 MHz s ' to -154 MHz s 



The selected groups of fibers (Table 1) were analyzed us- 
ing the wavelet method as described by Meszarosova et al. 
(2009a, b). We searched for tadpoles in the wavelet power spec- 
tra of radio flux time series and selected only dominant tadpoles 
corresponding to the 99% confidence level. The confidence level 
implies a test against a certain background level. If a peak in 
the wavelet power spectrum is significantly above background 
spectrum, then it can be assumed to be a true feature with a cer- 
tain percent of the confidence (Farge 1992). In the present pa- 
per the background spectrum is modeled by the red noise, and 
a computation of the confidence level is performed as described 
by Torrence & Compo (1998). The detected wavelet tadpoles 
and their parameters are summarized in Table 2. Their period P 
ranges from 1.1 to 2.3 s, and the drift of their head maxima FDt 
is in the range from -61 MHzs ' to -152 MHz s '. Comparing 
the frequency drifts of all fiber bursts in Table 2, we can see that 
the drift of fiber bursts is very similar to that of the head max- 
ima of the corresponding tadpoles (compare values in columns 3 
and 5 in Table 2). This can be also seen in Fig. [T] see the times 
tl = 12:01:13.7, r2= 12:01:13.1, f3 = 12:01:12.3, r4 = 8:35:15.5, 
f5 = 8:35:13.3, f6 = 8:35:10.8 UT at the selected frequencies 
1470, 1550, 1648, 1381, 1536, and 1733 MHz, respectively. 
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Table 1. Basic characteristics of the selected events with fiber bursts. 
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Nov. 23, 1998b 
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Nov. 23, 1998c 
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IN 
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Apr. 26, 2001b 
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Nov 18, 2003 


08:34 
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Jul. 11, 2005 
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Table 2. Parameters of the fiber bursts and associated wavelet tadpoles. 





Fiber bursts 


Wavelet tadpoles 


No. 


Frequency range 


FDp 


Period 


FDt 




[MHz] 


[MHz s-'] 


[s] 


[MHz s-'] 


1 


1456 ~ 2000 


-154 


1.5 


-152 


2 


1370- 1800 


-130 


1.5 


-122 


3 


1370- 1800 


-128 


1.4 


-127 


4 


1000- 1330 


-144 


1.1 


-136 


5 


1050- 1600 


-92 


2.0 


-87 


6 


1000- 1600 


-76 


2.3 


-61 


7 


1100-2000 


-73 


1.4 


-75 


8 


1100- 1800 


-78 


1.9 


-86 



Because the wavelet tadpoles are interpreted as signatures of 
the fast sausage magnetoacoustic wave trains (e.g. Nakariakov 
et al. 2004), the similarity of both the drifts indicates that the 
fiber bursts are physically connected with these magnetoacoustic 
wave trains. Based on these results we propose that the analyzed 
dm- fiber bursts are generated by the fast sausage magnetoacous- 
tic wave which modulates the radio emission of superthermal 
electrons trapped in a flare loop or in current sheet. 



3. Semi-empirical model 

This model is similar to the one presented by Kuznetsov (2006). 
However, it differs in one aspect: for the fast sausage magnetoa- 
coustic waves considered here, the magnetic and density wave 
perturbations {b, and p\) are in phase (Erdelyi, priv. communi- 
cation, 2012), contrary to Kuznetsov (2006), where these pertur- 
bations are in anti-phase. 

We assume that the magnetoacoustic wave is triggered by 
a flare and propagates upwards along a dense vertical flare loop. 
The plasma density and magnetic field in the loop are taken ac- 
cording to the models by Aschwanden (2002, 2004), where the 
magnetic field is expressed as B{h) - Boo (1 + hjhoY^ (h is 
the height in the solar atmosphere, ho = 75 Mm, and Boo is the 
footpoint magnetic field taken in the following runs No. 1 and 2 
as Boo - 65 G). The plasma density and magnetic field profiles 
used in the model are shown in Fig.|2] 

Although a real wave evolution is complex (see Roberts et al. 
(1984), Nakariakov et al. (2004), and also the foUowing MHD 
model), for a simplification we take the velocity of this wave 
perturbation as the Alfven speed va computed along the loop. 
The ratios of wave perturbations is taken as b^/Bo = pi/po 1, 
where Bo and po mean the local magnetic field and density in the 
Aschwanden's models. 



We took a wave in the form of the wave packet as 

\21 



F(h, f) - FqCo exp 



h - VA(h)t - hp 



X 



Co - exp 



xcos 



t - tAn 



2n(h - VA(h)t - ho) 



(1) 



dt 



(2) 



where F(h, t) means the density pi or magnetic field b^ wave per- 
turbation, F() is the wave amplitude, h the height in the loop, t the 
time, x'A the local Alfven speed, hp the wave initiation height, 
the spatial width of the wave packet, and di the wavelength. By 
changing ds for a fixed di we can generate the wave for a single 
fiber or group of fibers. We also added a factor Co that limits ra- 
dio emission to our chosen frequency interval 1300-1700 MHz 
(see the following): tAmax is thus the time of the maximum wave 
amplitude, and d, the corresponding characteristic time. In our 
computations the time t changes from to 5 s, fAmoA - 2.5 s, 
dt- 1.5 s, wdho- 16.9Mm.The other parameters are presented 
in Table 3. This propagating wave packet is then superimposed 
on the density and magnetic field profiles presented in Fig.|2] 

In agreement with Kuznetsov (2006), we assume that the ra- 
dio emission is produced at the double (harmonic) upper hybrid 
frequency as a result of a coalescence of two upper hybrid waves, 
generated by the loss-cone instability of superthermal electrons. 
Furthermore, we assume that the radio emission is generated 
with a constant intensity in the whole interval of considered 
heights. It means that every unit of the volume along the flar- 
ing loop (or current sheet) generates the same radio emission. 
Thus, the intensity of the radio emission at some specific fre- 
quency depends only on the gradients of and B in the radio 
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Fig, 2, Plasma density and magnetic field profiles in dependence 
on the height in the modeled loop. The horizontal full lines in 
the density plot mean the plasma densities corresponding to the 
maximum and minimum of frequencies considered in the radio 
spectra. The vertical dashed lines limit this density interval. 
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Fig. 4. Example of the time and frequency profiles of the fiber 
burst modeled in the semi-empirical model for run No. 1. The 
time profile is plotted at the 1500 MHz frequency and the fre- 
quency profile at the time 12.5 s; compare with the spectrum in 
Fig. 3. 
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Fig. 5. Computed radio spectrum with the fibers modeled in the 
semi-empirical model for run No. 2. 
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Fig. 3. Computed radio spectrum with the fiber modeled in the 
semi-empirical model for run No. 1 . 
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Fig. 6. Example of the time and frequency profiles of the fiber 
bursts modeled in the semi-empirical model for the run No. 2. 
The time profile is plotted at 1500 MHz frequency and the fre- 
quency profile at time 12.5 s; compare with the spectrum in 
Fig. 5. 



source; i.e., the wave perturbation modulates the radio emission 
as proposed by Treumann et al. (1990). Then, the resulting ra- 
dio flux at the specific frequency / can be expressed as a sum 
of emissions from plasma volumes having the upper hybrid fre- 
quency/up close to //2. 

First, we made two runs with the parameters shown in 
Table 3. The first run was performed for a single perturbation 
pulse and the second one for the wave train. The computed ra- 
dio spectrum is in Fig. [3] (Remark: For better presentation of all 
computed spectra at their beginning we added the 10-second ra- 
dio spectrum of the unperturbed state.) An example of the time 



Table 3. Computation parameters. 



Run No. 


Pi/Po 


b-JBo 


ds (km) 


dt (km) 


1 


0.01 


0.01 


200 


800 


2 


0.01 


0.01 


800 


800 



and frequency profiles of the resulting fiber for run No. 1 (cuts 
in the spectrum shown in Fig.[3]l is presented in Fig.|4] Besides 
the emission peak the emission depression on the low-frequency 
side of the fiber, as observed, is also clearly visible. The radio 
spectrum and time and frequency profiles for run No. 2, i.e. for 
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Fig. 7. Computed radio spectrum with the 10 fiber bursts mod- 
eled by the semi-empirical model. 



In our model we use the ideal 2D magnetohydrodynamic 
(MHD) equations: 



Dt ^ 
DB 



De 



Dt 



yeV ■ V, 



(3) 

(4) 

(5) 

(6) 
(7) 



the magnetoacoustic wave train, are shown in Figs. |5]and|6] As 
can be seen here, the wave train generates a series of fiber bursts. 

Furthermore, as an example, we simulated a group of ten 
fiber bursts. We consider the wave which is spatially limited to 
ten wavelengths (wavelength di = 800 km). It propagates up- 
wards in the solar atmosphere as in runs 1 and 2. However, 
the magnetic field Boo in the used magnetic-field model is 
Boo = 48 G. The constant Co (see Eq. 2) is taken as Co = 1, 
which means no global frequency modulation of the radio flux. 
The computed spectrum with ten fiber bursts is shown in Fig. Q 

Comparing this figure with the radio spectrum in the right- 
hand panel of Fig. [1] we can see some similarities as well as 
some deviations. Generally, using this semi-empirical model, 
we can simulate the frequency drift and repetition period of the 
fiber bursts by varying the Alfven speed (through the magnetic 
field Boo) and wavelength of the magnetoacoustic wave (c/^). 
On the other hand, deviations are caused by limited informa- 
tion about conditions in the real radio source and also by limita- 
tions in the semi-empirical model used. One of the main problem 
of this semi-empirical model is that dispersive properties of the 
magnetoacoustic waves are neglected. Therefore, in the follow- 
ing we describe a first attempt to simulate the fiber bursts in a full 
MHD model, where these properties are implicitly included. 



4. MHD model 

The main problem of any such model is an initiation of the suf- 
ficiently good waveguide for the fast sausage magnetoacoustic 
waves, which is in equilibrium in the gravitationally stratified 
solar atmosphere. For a waveguide we can use the dense loop 
or current sheet, see Jelmek & Karlicky (2012). In an analytical 
form, the model of the vertical current sheet in the gravitation- 
ally stratified atmosphere proposed by Galsgaard and Roussev 
(2002) is an appropriate one for initiating the equiUbrium waveg- 
uide. 

In our case the vertical current sheet (it would be better to 
call it the current layer) is only a prototype of the waveguide 
guiding the magnetoacoustic waves. We did not consider any re- 
connection and fragmentation processes in this current sheet as 
discussed in Barta et al. (2011a, b), Karlicky et al. (2012), and 
Cargill et al. (2012), among others. Although the fragmentation 
is real in these current sheets, we think that such vertical cur- 
rent sheets can exist; see, e.g., the observed current sheet and its 
width on the order of 1000 km in Fig. 3 in the paper by Lin et al. 
(2005). 



where D/Df = d/dt + v-V is the total time derivative, p is a mass 
density, v flow velocity, B the magnetic field, and g = [0, -gQ, 0] 
is the gravitational acceleration with gQ - 214 m ■ s"^. The 
current density j in Eq. (|4]i is expressed as 



J = — V X B, 



(8) 



where //o = 1.26 x 10 ^ H ■ m is the magnetic permeability of 
free space. The specific internal energy e in Eq. (|6]l is given by 

(9) 



(r-i)e' 

with the adiabatic coefficient y - 5/3. 

Equations ([3]l-(|6|l are solved numerically by the FLASH code 
(Fryxell et al. 2000, Lee & Deane 2009). This code implements 
second- and third-order unsplit Godunov solvers and adaptive 
mesh refinement (AMR) (see e.g. Chung 2002, Murawski 2002). 
Spatial resolution of the numerical grid is determined by the 
AMR method, and we use the AMR grid with a minimum (max- 
imum) level of refinement blocks set to 3 (6). The total num- 
ber of computational cells covering a whole simulation plane 
is 267456. A spatial cell size has to be much smaller than the 
width of the current sheet in the x-direction and the minimal 
wavelength of the magnetoacoustic waves in the h-direction. 
In our model the minimal cell sizes are Ax = 0.0175 Mm and 
Ah - 0.005 Mm (compared with the half width of the current 
sheet wcs = 0.35 Mm and the minimal wavelength of about 
1.0 Mm). 

At the start of our numerical calculations we set a refinement 
procedure in the region covering the current sheet. During the 
calculations, the FLASH code automatically controlled in each 
time step the gradient of mass density, which value is then used 
for refining the grids. For our numerical simulations we used 
the two-dimensional (2-D) Eulerian box (-2.25,2.25) Mm x 
(13.0, 23.0) Mm in the x- and h- (height) directions, respectively. 
As a consequence of an extension of the real plasma medium we 
apply free-boundary conditions at the boundaries of the simula- 
tion region, so that the waves can freely leave the simulation box 
without any significant reflection. 

For the description of magnetic field in the vertical current 
sheet in gravitationally stratified solar atmosphere in the x - h 
plane we use following expressions: 



B J x,h)^B,^ In 
A 



cosh I — 

Wcs 



h 

expl-- 



Bhi^x, h) - Bo tanh I — | exp ( - - 



(10) 



(11) 
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where is a magnetic field at x — > oo, wcs is the half width of 
the current sheet, and A denotes the magnetic scale-height. For 
an equilibrium (v = 0), the Lorentz and gravity forces have to be 
balanced by the pressure gradient in the entire physical domain 



-Vp+jxB + gg ^0. 



(12) 



From this condition we can derive the formulae for distribu- 
tion of the mass density and gas pressure. For the details see 
Galsgaard & Roussev (2002) and Jelmek et al. (2012). 

In the initial state, we generate the vertical current sheet with 
the density profile along its axis according to the Aschwanden 
model (2002). The half width of the current sheet at the height 
h = 16.9 Mm is taken as wcs - 0.35 Mm. In the paper by Jelmek 
& Karlicky (2012), we found that the neutral current sheet serves 
as the same waveguide as the dense slab if the Alfven speed at 
the half width of the current sheet is the same as in the slab. 
Therefore, in the present current sheet we took this Alfven speed 
as 690 km s ' . The change of this speed in the studied interval of 
heights in the current sheet is small. 

At the start of the numerical simulation, the initial equilib- 
rium state is perturbed by the Gaussian pulse in the x-component 
of velocity and has the following form (e.g. Nakariakov et al. 
2004, 2005): 



(13) 



X 


(x-/2p)2 






— exp 


[ \ 


■ exp 





where Aq is the initial amplitude of the pulse, and A^- - Ah — 
0.35 Mm are the widths of the velocity pulse in the longitudinal 
and transverse directions. This pulse tends to trigger fast sausage 
magnetoacoustic waves. The initial perturbation is located at the 
height /ip - 15 Mm. 

Using this MHD model, we computed an evolution of the 
density profiles along the axis of the vertical current sheet (see 
examples in Fig. |8]l. These density profiles were then used in 
computations of the artificial radio spectrum, as well as the time 
and frequency profiles shown in Figs. |9] and [10] We used the 
same simple radio emission model as in the case of the semi- 
empirical model. As seen in Fig. |9]the frequency drift of these 
fiber bursts partly changes. It is caused by the wave dispersive 
effects (included in this MHD model), which change the density 
wave profile during the wave propagation. Compare this with 
the semi-empirical model, where the perturbation profile is rigid. 
The maximum of fiber bursts are produced at the locations with 
minimal density gradients and these gradients evolve in time. 
Similarly, the emission depressions, corresponding to the maxi- 
mal density gradients, evolve in time. 

The detailed analysis of these results shows that the model 
used for the vertical current-sheet is not the best waveguide for 
the studied magnetoacoustic waves. The change in the mag- 
netic field gradient across the current sheet is not sharp enough. 
Therefore, some wave energy escapes from this waveguide. This 
wave energy leakage limits the wave train length and thus also 
the number and duration of computed fibers. We hope that in 
a future model we will generate a better waveguide, in which 
a longer wave train and thus more fiber bursts will be generated 
by a single perturbation. On the other hand, more fiber bursts can 
be generated by a time repetition of the initial perturbation be- 
cause this is also possible in real conditions. Furthermore, some 
real radio spectra might show a superposition of the fiber bursts 
generated in several nearby waveguides. 




5x10=" 
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h (km) 

Fig. 8. The density profiles in the MHD model at 1 1 s (full line) 
and 14 s (dash-dotted line). The times correspond to times in the 
spectrum in Fig. 9. 
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Fig. 9. The computed radio spectrum with the fibers modeled in 
the MHD model. 
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Fig. 10. An example of the time and frequency profiles of the 
fiber bursts modeled in the MHD model. The time profile is plot- 
ted at the 1500 MHz frequency and the frequency profile at the 
time 14.5 s; compare with the spectrum in Fig. 9. 



5. Discussion and conclusions 

After analyzing eight groups of observed fiber bursts by the 
wavelet technique, we found a connection between the dm- 
fiber bursts and fast sausage magnetoacoustic waves. Based on 
this connection, we developed the model of the dm-fiber bursts 
with the fast sausage magnetoacoustic wave propagating in the 
waveguide of the enhanced plasma density: (a) in the dense loop 
according to the Aschwanden's density and magnetic field mod- 
els (semi-empirical model) and (b) in the vertical and gravita- 
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tionally stratified current-sheet (MHD) model. In both of these 
models, using the simple emission model according to Treumann 
et al. (1990), we computed the artificial radio spectra with the 
fiber bursts that resembles to the observed ones. 

We assumed that every unit of the volume along the flar- 
ing loop (or current sheet) generates the same radio emis- 
sion. The intensity of the radio emission at some specific fre- 
quency depends only on the density gradient in the radio source. 
However, the wave perturbation (with weak local magnetic mir- 
rors) can also modulate the growth rate of the upper hybrid 
waves through change in the loss-cone angle of superthermal 
electrons (Winglee & Dulk 1986, Yasnov & KarUcky 2004) and 
thus change local emissions. For this purpose, the weighting 
function Weih{t)) that expresses these variations can be included 
in future models. 

In all studied decimetiic fiber bursts, the frequency drift of 
the wavelet tadpoles coiTesponds to the drift of individual fiber 
bursts. It speaks in favor of the present model with the magnetoa- 
coustic waves in comparison to the model with whistler waves, 
which is usually considered in the metric frequency range. We 
propose to do similar wavelet analysis of the metric fiber bursts 
and to check the validity of their models. 

While in the semi-empirical model the form of the wave per- 
turbation is rigid, the MHD model describes perturbations in- 
cluding dispersive properties of the magnetoacoustic waves. It 
leads to a wavy character of the fiber bursts, which can even 
be seen in some observed fiber bursts (see e.g. the fibers in the 
upper-left part of Fig. [!}. 

The semi-empiiical model is very similar to that of 
Kuznetsov (2006), except that the density and magnetic field 
perturbations are in phase. The phase relation presented in 
Kuznetsov (2006) is not correct for the fast sausage magne- 
toacoustic waves. However, this difference in both the mod- 
els leads to only small differences in the radio spectrum, be- 
cause the electron-cyclotron frequency coce is much lower than 
the plasma frequency ojpe, and thus the upper-hybrid frequency 
roughly equals the plasma frequency and the magnetic field can 
be neglected in this case. However, the variation of the magnetic 
field in the magnetoacoustic wave train can influence the growth 
rate of the upper hybrid waves owing to the loss-cone instabil- 
ity of superthermal electrons. It can modify the radio emission 
at some locations along the wave train and thus change the fiber 
burst profiles. This effect is not considered in the present paper. 

Comparing both these models, the semi-empirical is easy to 
use, however, with a limited description of physical processes 
under study. The MHD model is much more complicated. One 
of the main problems is to generate the appropriate initial equi- 
librium state. The MHD model includes the dispersive effects of 
the propagating magnetoacoustic waves, so it describes the fiber 
burst generation in a more realistic way. However, the present 
MHD vertical and gravitationally stratified waveguide (the cur- 
rent sheet according to the Galsgaard and Roussev model, 2002) 
is not ideal. Owing to its relatively smooth boundaries, the en- 
ergy of the magnetoacoustic waves escape from this wave guide. 
The energy of wave decreases, and it is therefore difficult to 
make long-lasting wave train. This means that we can only sim- 
ulate a few fiber bursts per one perturbation instead of the long 
series as in the semi-empirical model. This problem is connected 
to a solution of the initial equilibrium waveguide. Therefore, to 
make this MHD model of the fiber bursts more realistic, a better 
initial waveguide needs to be found. On the other hand, the series 
of fiber bursts can be generated in the present MHD model us- 
ing time series of perturbations. Furthermore, some of real radio 



spectra might show a superposition of the fiber bursts generated 
in several nearby waveguides. 

Considering our simple emission model, we plan to solve 
the inverse problem to the presented solution: i.e., we plan to 
determine the density profiles of the propagating magnetoacous- 
tic wave from the fiber burst profiles measured along the radio 
frequency at some specific times. 
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